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Hexagonal HoMnOa is a frustrated antiferromagnet (Tjv=72 K) ferroelectric (Tc7=875 K) in which 
these two order parameters are coupled. Our neutron measurements of the spin wave dispersion for 
the = 2 Mn^+ on the layered triangular lattice are well described by a two-dimensional nearest- 
neighbor Heisenberg exchange J=2.44 meV, and an anisotropy D that is 0.093 meV above the 
spin reorientation transition at 40 K, and 0.126 meV below. For H \\ c the magnetic structures 
and phase diagram have been determined, and reveal additional transitions below 8 K where the 
ferroelectrically displaced Ho^^ ions are ordered magnetically. 

PACS numbers: 75.30.Ds, 75.30.Et, 75.50.Ee 
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Magnetic ferroelectric materials are rare, and rarer still 
are systems where these two disparate order parameters 
exhibit significant coupling 0. Such multiferroics have 
been of particular interest recently both to understand 
the fundamental aspects of this coupling 2, 3, 
and because of the intriguing possibility of using these 
coupled order parameters in novel device applications 
• The hexagonal HoMnOs system of particular interest 
here is a prototype multiferroic where the Ho-0 displace- 
ments give rise to a ferroelectric moment (Tq = 875K) 
along the c-axis, the Mn moments order at 72 K, and the 
order parameters are naturally coupled through the Ho- 
Mn exchange and anisotropy interactions. The magnetic 
system has the added interest that the Mn moments oc- 
cupy a fully frustrated triangular lattice. Our inelastic 
neutron scattering measurements of the Mn spin dynam- 
ics reveal the planar nature of the spin system, and allow 
us to establish the basic model for the magnetic inter- 
actions in the system and determine the interaction pa- 
rameters. Field-dependent neutron diffraction measure- 
ments reveal the nature of the magnetic phase diagram, 
and in particular the Ho involvement in the hysteretic 
magnetic transitions at low temperatures. Our results 
also suggest that the Ho coupling may be responsible for 
the first-order spin-flop transition found at intermediate 
temperatures. 

We have grown single crystals of HoMnOa using a 
traveling solvent optical floating zone furnace, and the 
diffraction and inelastic neutron measurements were per- 
formed at the NIST Center for Neutron Research on the 
BT2 and BT9 thermal triple-axis instruments. The mag- 
netic structure of HoMnOa has been studied previously 
with neutron powder diffraction |lCl . optical second- 
harmonic generation (SHG) |iull2.ll3 . ll4l |. and magnetic 
and dielectric susceptibility The Mn^+ ions are ar- 
ranged in a 2-dimensional (2D) triangular lattice, with 
successive layers offset from each other. Fig. shows 
the integrated intensities of the (1,0,0) and (1,0,1) (com- 
mensurate) magnetic Bragg peaks, where we see that be- 
low Tn « 72K the spin-2 Mn'^^ moments order anti- 
ferromagnetically in a non-coUinear 120-degree structure 
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FIG. 1: (color online). Neutron diffraction measurements of 
the (1,0,0) and (1,0,1) magnetic Bragg reflection integrated 
intensity. Two spin reorientation transitions (indicated by 
dashed lines) lead to changes in the intensities of both peaks. 
The Mn^+ spin configurations for these phases are shown in 
the schematics to the right, where open circles indicate the 
position of Mn ions ai z — 0, filled circles indicate Mn ions 
a.t z — c/2, and arrows indicate the direction of the local 
magnetic moment. Ho'^"'' moments (not shown) order in the 
low-temperature phase. 



within each plane. At Tsr ~ 40 K a sharp spin reorien- 
tation transition is observed associated with a change in 
the magnetic symmetry to the Pd'^cm' structure, as in- 
dicated in Fig. n] Below 8 K half of the Ho^+ moments 
order [9lll0l|. which results in a third zero-field spin reori- 
entation transition of the Mn spins. The three zero-fleld 
Mn spin structures are also shown in Fig. ^ 

Measurements of the spin-wave dispersion were per- 
formed on a 1.8 gram single crystal in both the (H,K,0) 
and (H,0,L) scattering planes. An example of a constant- 
Q scan at 20K showing two distinct magnetic modes is 
shown in Fig. Fits to individual scans were used to 
establish the peak positions. The dispersion relations at 
20K, in the intermediate-temperature phase, are plotted 
in Fig. I^b and [2t along two different directions in re- 
ciprocal space. Data for the spin-wave dispersion along 
the (1,0, L) direction (inset Fig. ^) show no discernible 
dispersion in the out-of-plane direction, indicating pri- 
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FIG. 2: (color online). Spin-wave dispersion at 20K. (a) Mea- 
surements of the spin waves at Q=(1.3,0,0). Data were taken 
on BT2 with horizontal collimations of 60'-40'-40'-80' and a 
fixed final energy of 14.7 meV. The line shows a fit to the 
data using two modes convoluted with the resolution func- 
tion. Peak positions from both constant-energy and constant- 
Q scans are plotted for two in-plane directions in reciprocal 
space in (b) and (c), along with the three modes of the dis- 
persion curve (solid lines) given by Eq. |1]( J = 2.44(7) meV 
and D = 0.126(25) meV). Dashed lines indicate two (disper- 
sionless) crystal field levels of Ho at 1.5(1) and 3.1(1) meV. 
Inset: L dependence of the spin-wave spectrum, showing no 
discernable dispersion. 



marily 2D spin dynamics. Data were also taken at 50K, 
in the high-temperature phase, and a comparison of two 
constant-Q scans at the magnetic zone center are shown 
in Fig. 13^. The zone-center gap decreases significantly, as 
shown in the comparison of the 20K and 50K dispersion 
data in Fig. ^p. 

The dispersion of the excitations in HoMnOa was cal- 
culated using a linear spin-wave analysis. We start from 
the Hamiltonian 



(1) 



where indicates the sum is over nearest- neighbor 

in-plane pairs, J is the primary antiferromagnetic ex- 
change, and D is the anisotropy. Introducing three dif- 
ferent flavors of bosons and using a Holstein-Primakoff 
transformation, we derive a linearized Hamiltonian which 
we diagonalize by mapping to the Hamiltonian of the 
quadratic quantum mechanical oscillator, similarly to Ref 
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FIG. 3: (color online), (a) Inelastic neutron measurements 
at Q=(2,0,0) at 20K and 50K. Data were taken on BT9 with 
horizontal collimations of 40'-40'-40'-40' and fixed incident 
energy of 13.7 meV. Lines show fits to the data. The increase 
in intensity at lower energies is primarily due to a crystal field 
excitation at 3.1 meV. (b) Comparison of 20K and 50K dis- 
persion data near the zone center. Data have been reduced to 
the first Brillouin zone. Dotted lines show spin-wave disper- 
sion curves for 20K data, as in Fig. |21 Solid lines show fit of 
Eq. Blto 50K data (J = 2.44 meV is fixed and D = 0.093(10) 
meV). 



Il5l. We define the lattice Fourier sum as 
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For simplicity, we now define 



Zl = z + z 

Z2 = -{z + z*)/2 + i{V3/2){z- z*) 

Z3 ^ -{z + z*)/2-i{V3/2){z~z*). (3) 



We obtain three modes whose dispersion is given by 



= 35Jv/(l + + 2D/J) (1 - 2z,), 



(4) 



where i = 1, 2, or 3, and S* = 2 is the spin at the Mn'^^ 
ion. Equation^was used to fit our data, with the results 
shown by the curves in Figs. [2b; Eti and|2l3. This sim- 
ple model provides a remarkably good description of the 
spin dynamics. Quantitatively, we obtain J = 2.44(7) 
meV and D = 0.126(25) at 20 K. Data taken away from 
the zone center showed little temperature dependence, 
while at the zone center there was a significant difference 
between the 20K and 50K data, indicating that the tem- 
perature dependence of the spin-wave spectrum comes 
largely from changes in the anisotropy. For the 50K 
data, J was therefore kept fixed at 2.44 meV, yielding 
D — 0.093(10). The increase in the anisotropy D at 
lower temperatures is signficant, and suggests that this 
anisotropy, likely originating from the holmium, drives 
the spin reorientation transition. 
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FIG. 4: (color online), (a) (1,0,0) magnetic Bragg peak inten- 
sity versus magnetic field for different temperatures (offset by 
15,000 counts for clarity), (b) (2,1,0) magnetic Bragg peak 
intensity versus magnetic field for different temperatures. Cir- 
cles are data taken with increasing field, while triangles are 
data taken with decreasing field. Dashed lines show 0.4K 
transitions, with two intermediate phases a and c and hys- 
teretic overlap in region b. 



In the presence of a magnetic field Tsb. shifts to lower 
T and broadens, and HoMnOa develops a reentrant 
PG'^c'm phase. A sharp anomaly in the dielectric sus- 
ceptibility at TsR indicates a coupling between the mag- 
netic and ferroelectric order Q . SHG measurements as a 
function of magnetic field also show the change in mag- 
netic symmetry and the reentrant phase as a function of 
magnetic field suggested by dielectric susceptibility [l3l |. 
More recent measurements show a marked decrease in 
the strength of the SHG signal (which is due to the mag- 
netic ordering of the Mn'^+ ions) when an electric field is 
applied 0. 

To investigate the magnetic field-dependent phase di- 
agram, a 0.5 gram sample was mounted in the (H,K,0) 
scattering plane inside a vertical-field 7 Tesla supercon- 
ducting magnet with a Helium-3 insert. A magnetic 
field of up to 5 Tesla was applied along the c axis, and 
the (commensurate) (1,0,0) and (2,1,0) magnetic peaks 
were used to track the evolution of the magnetic phases 
with magnetic field and temperature. In the interme- 
diate temperature phase, a sufficiently strong applied 
magnetic field along the c axis pushes HoMnOa into the 
high-temperature phase 0, 0| . The boundary between 
the intermediate- and high-temperature phases (shown in 
Fig. [SI as a function of temperature and magnetic field) 
is in agreement with the phase diagram established from 
dielectric susceptibility measurements Il6j . 
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FIG. 5: (color online). Temperature and magnetic field phase 
diagram for HoMnOa obtained from diffraction measurements 
of the (1,0,0) and (2,1,0) Bragg reflections. Lines are guides 
to the eye. Regions a and c are intermediate phases, with 
hysteretic overlap in region b, as in Fig. 0] Region d is an in- 
termediate region where peak intensities change continuously. 
Dashed line is approximate location of low-temperature Ho'^"'" 
ordering transition, and dotted line is approximate transi- 
tion observed by SHG and magnetic susceptibility measure- 
ments [73 . [T^ . These transitions were not observable from 
the (1,0,0) and (2,1,0) diffraction data. 



At lower temperatures the transition becomes consid- 
erably more complicated. Figure 0Jl shows the peak in- 
tensity of the (1,0,0) Bragg scattering as a function of 
magnetic field at several temperatures. At higher tem- 
peratures, the intensity decreases approximately linearly 
with increasing field over a well-defined transition re- 
gion. At 0.4K, however, there are two distinct, step-like 
transitions, with considerable hysteresis in the lower-field 
transition, indicating the existence of a well-defined in- 
termediate phase between the low- and high-field phases. 
As the temperature increases, these transitions become 
broader and the hysteresis decreases, merging together 
into the single transition (without appreciable hystere- 
sis) observable at higher temperatures. 

Figure^t shows equivalent measurements of the (2,1,0) 
Bragg peak intensity. Again, at high temperatures, an 
approximately linear change in scattering intensity over 
a well-defined transition region is observable. At low tem- 
peratures, as for the (1,0,0) refiection, the situation be- 
comes more complicated. Instead of two, we find three 
transitions, with an additional lower-field transition not 
observable in the (1,0,0) data. Unlike the (1,0,0) refiec- 
tion, the scattering intensity of the (2,1,0) refiection does 
not change monotonically. 

Figure O shows a phase diagram extracted from the 
(1,0,0) and (2,1,0) diffraction data. Below 4K there 
are two distinct intermediate phases (labeled a and b) 
between the zero-field and high-field phases, with signifi- 
cant hysteretic overlap between them. These data are not 
sensitive to the low-temperature spin reorientation ob- 
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servable in the (2,1,0) difFraction data (Fig. indicated 
roughly by the dashed line in Fig. |S1 It is interesting 
to note, however, that the intermediate phases are only 
observable below the low-temperature spin reorientation 
transition where Ho^+ ordering becomes significant. 

At higher temperatures, our phase diagram agrees well 
with previous dielectric measurements '3]. These SHG 
measurements did not observe any intermediate phase at 
low temperatures while more recent magnetic sus- 
ceptibility and heat capacity measurements have iden- 
tified one intermediate phase While we observe 
two unambiguous, distinct intermediate phases at low 
temperatures, the combined intermediate phases we ob- 
serve agree qualitatively with the intermediate phase re- 
cently observed Further measurements of out-of- 
plane magnetic peaks will be necessary in order to deter- 
mine the exact magnetic structure of these intermediate 
phases. Recent SHG measurements suggest that Ho'^"'" 
magnetic moments can even be ordered ferromagnetically 
with the application of an electric field . Neutron scat- 
tering reveals considerable complexity in the magnetic 
field phase diagram for HoMnOs, and extending these 
neutron measurements to the electric field phase diagram 
may also prove similarly revealing. 

The out-of-plane magnetic interactions in HoMnOs, 
both between successive Mn'^+ layers and between Mn'^+ 
and Ho'^+, are almost completely frustrated. It is there- 
fore not surprising that the spin dynamics in the interme- 
diate and high temperature phases are two-dimensional 
in nature, and we find that the spin- wave spectrum is well 
described by Eq. Ferroelectric lattice distortions re- 
lieve the interplane frustrations in HoMnOa '3| , and may 
be the driving force behind the spin-reorientation tran- 



sitions of the Mn^+ magnetic lattice. The role of Ho in 
these reorientation transitions may be very similar to the 
role of Nd in Nd2Cu04, another layered transition- metal 
oxide in which frustrated interplanar interactions enable 
weak, higher-order interactions to dictate the magnetic 
structure. Nd2Cu04 has a non-coUinear structure ^3 
with Cu spin-reorientation transitions [isl ll^ that are 
controlled by the rare earth single- ion anisotropy pol[2ll |. 
We suggest that it is a similar interaction in HoMnOa 
that controls the spin-flop transition, and is the origin of 
the magnetic-ferroelectric coupling. 

There has been considerable experimental progress in 
recent years on mapping out the magnetic phase diagram 
of HoMnOa and discovering signatures of the interplay of 
magnetic and ferroelectric order in various physical prop- 
erties. However, a theoretical understanding of these 
interactions is still lacking. Although a complete the- 
oretical description of HoMnOs will necessarily require 
significantly more complexity than Eq. ^ having exper- 
imentally determined the dispersion relations and values 
for the primary magnetic exchange constants is an im- 
portant step toward being able to model the interaction 
between ferroelectric and magnetic order in HoMnOa. 
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